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ABSTRACT 

The  problem  of  estimating  the  acoustic  field  produced  by  a point  source 
located  in  an  environment  which  is  arbitrarily  stratified  with  depth  only  has 
received  considerable  attention.  Various  models  exist  which  provide  reasonable 
results  provided  sufficient  information  about  the  environment  and  source  char- 
acteristic is  available.  This  paper  is  concerned  with  solving  for  the  inverse 
of  this  problem,  i.e.,  information  about  the  environment,  in  particular  the 
impedance  at  the  ocean-bottom  interface  is  sought  from  a knowledge  of  the  re- 
ceived acoustic  field  and  the  source  which  produced  it.  The  theoretical  for- 
mulation of  the  scheme  is  presented  and  its  numerical  implementation  is  dis- 
cussed and  illustrated  for  three  examples  of  successive  orders  of  environmental 
complexity. 
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Although  the  subject  of  this  paper  is  the  inverse  solution  we  will  begin 
with  a discussion  of  the  direct  solution  because  it  came  first  chronologically 
and  its  success  provided  the  motivation  for  the  development  of  the  inverse 
solution. 

By  direct  solution,  figure  1,  we  mean  that  situation  where  the  source 
and  environmental  characteristics  are  assumed  to  be  known  and  a solution  for 
the  received  field  is  desired.  The  two  major  assumptions  for  both  the  direct 
and  inverse  solutions  are:  that  the  source  is  concentrated  at  a single  point 
in  space  in  an  ocean  with  cylindrical  symmetry,  and  that  the  sound  speed  is 
allowed  to  vary  arbitrarily  with  depth  but  is  invariant  with  range. 

It  is  well  known  that  the  field  can  be  obtained  by  performing  a Fourier 
synthesis  over  frequency  of  the  product  of  H(r,z,f);  the  transfer  function 
of  the  medium,  and  F(f),  which  is  the  spectrum  of  the  source  waveform.  It  is 
also  well  known  that  the  transfer  function  can  be  expressed  as  the  Fourier- 
Bessel  transform  given  by  the  bottom  equation  where  G is  the  depth  dependent 
Green's  function  which  contains  all  of  the  information  about  the  environment. 

Both  of  the  integrals  can  be  numerically  evaluated  using  the  discrete 
fast  Fourier  transform,  commonly  known  as  the  FFT.  The  evaluation  of  the  first 
equation  is  rather  straightforward,  but  the  second  requires  some  explanation. 

The  transfer  function  given  in  figure  2 is  first  rewritten  in  terms  of 
the  Hankel  function  corresponding  to  outward  propagating  waves.  The  key  obser- 
vation is  that  if  the  large  argument  asymptotic  approximation,  given  by  the 
second  equation,  is  substituted  for  the  Hankel  function  the  transfer 

function  is  then  given  by  the  Fourier  type  integral  shown  on  the  third  line. 

It  should  be  noted  that  the  range  of  integration  is  relatively  quite  small 
since  it  is  related  to  the  minimum  sound  speed  of  the  problem.  We  have  in- 
cluded the  parameter  e to  indicate  that  in  practice  we  integrate  to  a slightly 
larger  wave  number  to  insure  that  evanescent  waves  are  also  included. 

The  third  equation  is  then  numerically  evaluated  using  the  FFT  for  each 
discrete  frequency  of  interest  to  provide  the  transfer  function  of  the  medium 
as  a function  of  both  range  and  frequency  at  a given  source  and  receiver  depth. 
Multiplication  by  the  source  spectrum  and  the  evaluation  of  another  FFT  then 
yields  the  received  waveform  as  a function  of  range  and  time. 

The  success  experienced  with  this  approach  led  us  to  consider  the  possi- 
bility of  inverting  the  entire  process. 

We  define  the  inverse  solution  given  in  figure  3,  as  that  situation  where 
the  source  characteristics  and  received  signal,  as  a function  of  range,  are 
assumed  to  be  known  and  information  about  the  environment  is  desired.  In  par- 
ticular, the  information  sought  is  the  Impedance  at  the  water  bottom  interface 
as  a function  of  wave  number. 

The  inverse  solution  is  obtained  by  first  performing  a deconvolution 
(i.e.,  taking  the  ratio  of  the  Fourier  transform  of  the  output  to  the  input) 
and  then  taking  the  inverse  Bessel  transform  over  horizontal  range  to  give  us 
the  depth  dependent  Green's  function  versus  the  horizontal  component  of  the 
wave  number.  Since  this  function  depends  on  the  properties  of  the  medium,  it 
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is  possible  to  carry  the  process  one  step  further  and  solve  for  the  impedance 
at  the  water  bottom  interface  under  the  assumption  that  you  know  the  sound 
speed  to  that  depth. 

Before  proceeding  to  the  examples  let's  examine  the  inversion  of  the 
transfer  function  itself  shown  by  the  first  equation.  In  figure  4.  As  before 
use  is  made  of  the  asymptotic  expression  for  the  Hankel  function  to  yield  the 
Fourier  type  integral  on  the  second  line.  The  effect  of  this  approximation  is 
minimal  and  will  be  illustrated  in  the  examples.  As  can  be  seen  from  the  last 
equation,  the  approximate  Green's  function  is  given  by  the  convolution  of  the 
exact  Green's  function  with  the  Fourier  transform  of  the  range  window  which 
extends  from  0 to 

Thus,  if  the  inverse  result  is  to  be  a good  approximation  to  the  exact 
answer  over  all  values  of  the  horizontal  wave  number,  the  aperture  in  range 
or  Rmax  must  be  fairly  large.  This  is  in  contrast  to  the  small  wave  number 
aperture  Cmax  needed  in  the  direct  solution.  The  effects  of  various  range 
apertures  will  be  graphically  presented  in  the  last  example  to  be  considered. 

It  is  interesting  to  note  that  the  approximate  inverse  solution  as  given 
by  the  middle  equation  can  be  considered  as  the  output  of  a long  horizontal 
line  array  with  a phase  delay  equal  to  the  horizontal  component  of  the  wave 
number  and  amplitude  shading  equal  to  the  square  root  of  the  distance  from  the 
source. 

The  results  of  inverse  solution  for  three  problems,  figure  5,  of  successive 
orders  of  environmental  complexity  will  be  presented.  In  all  cases  the  direct 
solution  was  used  to  generate  the  received  pressure  field.  Thus,  the  exact 
answer  for  the  inverse  solution  is  known  in  advance  and  these  results  serve  as 
a check  on  the  computations  and  as  an  illustration  of  the  effects  of  using  a 
limited  range  aperture. 

The  first  case  is  the  simplest  possible  example.  The  ocean  is  infinite 
and  of  constant  sound  speed  with  a CW  signal  as  the  input  waveform.  The  direct 
path  will  be  inverted  to  obtain  the  total  depth  dependent  Green's  function. 

The  next  case  is  slightly  more  complicated  since  a high  speed  half  space 
has  been  inserted  and  the  input  waveform  is  a pulsed  CW.  The  portion  of  the 
field  containing  the  reflected  and  lateral  waves  will  be  inverted  to  first 
yield  the  total  Green's  function  and  from  that  the  solution  for  the  Rayleigh 
plane  wave  reflection  coefficient. 

The  last  example  is  the  most  complicated  environment  that  can  be  considered 
under  the  initial  assumption  of  range  invariance.  The  sound  speed  profile  was 
taken  from  a deep  water  site  in  the  North  Atlantic,  and  numerous  multipaths 
are  present  in  the  received  signal.  The  received  pressure  field  for  an  infinite 
CW  input  waveform  will  be  inverted  to  produce  the  total  depth  dependent  Green's 
function. 

The  exact  inverse  solution  for  the  case  of  only  a direct  path  is  given  by 
the  first  equation  in  figure  6.  It  should  be  noted  that  it  is  necessary  to 
include  attenuation,  represented  by  the  parameter  d in  the  received  field,  in 
order  to  achieve  a practical  range  aperture. 
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The  Bessel  transform  on  the  right  was  evaluated  by  using  the  FFT  technique  , 

discussed  earlier.  The  associated  parameters  are  shown  on  the  bottom  of  the 
figure.  The  calculation  yields  values  of  the  Green's  function  for  wave  numbers 
from  to  to  ^max  w^^ch  Is  larger  than  the  free  space  wave  number  ki. 

The  Inverse  solution  for  the  depth  dependent  Green's  function  and  the 
exact  answer  are  shown  In  figure  7,  plotted  versus  the  horizontal  component  of  * 

the  wave  number.  The  agreement  Is  excellent  (the  two  curves  over  plot  exactly 
except  for  values  of  the  horizontal  component  of  the  wavenumber  less  than  .01. 

The  heavy  line  Is  the  IFFP  solution.)  except  for  very  striall  values  of  the  wave 

number  which  correspond  to  small  angles  as  measured  from  the  vertical.  For 

these  angles  the  horizontal  range  would  also  be  small  and  In  the  limiting  case  i 

of  the  receiver  directly  over  the  source  both  would  be  zero.  Thus,  the  product 

of  the  two  quantities  Is  also  small  and  violates  the  large  argument  assumption 

used  to  approximate  the  Hankel  function.  This  limitation  can  be  avoided  by  I 

evaluating  the  Bessel  transform  with  a FFT  but  In  a slightly  different  manner. 

This  has  been  done,  but  the  lack  of  time  does  not  permit  a more  detailed  ex- 
planation. 

The  received  waveforms  for  the  second  example  as  calculated  by  the  direct 
solution  for  only  the  bottom  Interfacing  portion  of  the  field  are  shown  In 
figure  8 for  three  ranges.  From  left  to  right  they  correspond  to  a range  prior 
to  the  critical  range  where  the  grazing  angle  Is  39. 60,  and  Is  greater  than 
the  critical  angle.  At  this  range  only  the  reflected  wave  exist;  next  for  a 
range  just  beyond  the  critical  range  where  the  grazing  angle  Is  31.7°  and 
slightly  smaller  than  the  critical  angle.  At  this  range  both  the  reflected 
and  lateral  waves  exist  but  are  coincident  In  time;  and  finally  at  a range  well 
beyond  the  critical  range  where  the  grazing  angle  Is  6.8°  and  the  two  waves 
are  separated  In  time.  These  signals  and  those  at  other  ranges  were  then  used 
as  Input  values  for  the  inverse  solution. 

The  right-hand  side  of  the  equation  In  figure  9 Is  evaluated  numerically 
by  first  taking  the  FFT  over  time  of  the  received  waveforms  at  each  range  point 
which  gives  the  function  script  G under  the  Integral.  The  Bessel  transform  Is 
then  evaluated  by  the  FFT  over  range  to  provide  values  for  the  depth  dependent 
Green's  function  appearing  on  the  left-hand  side  of  the  equation  versus  wave 
number.  Since  the  analytic  expression  for  G,  given  by  the  second  equation.  Is 
known,  we  can  then  solve  for  the  Rayleigh  plane  wave  reflection  coefficient, 

V(C). 

Before  looking  at  those  results.  It  should  be  noted  that  the  maximum  range 
of  Integration  used  In  this  example  corresponds  to  an  Incident  angle  of  83° 
so  that  the  range  aperture  Is  truncated  and  some  aliasing  In  the  final  result 
Is  expected. 

The  absolute  magnitude  of  the  Inverse  solution  and  the  exact  answer  are 
plotted  versus  the  angle  of  Incidence  In  figure  10.  It  Is  evident  that  excel- 
lent agreement  exist  everywhere  except  for  small  angles  less  than  2°,  and  for 
large  angles  greater  than  75°.  Small  angles  of  Incidence  are  associated  with 
small  values  of  horizontal  range,  thus  the  product  of  the  horizontal  component 
of  the  wave  number  and  range  Is  also  small;  and  this  violates  the  large  argu- 
ment assumption  used  In  approximating  the  Hankel  functions. 
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The  oscillation  at  the  large  angles  of  incidence  is  due  to  the  aliasing 
introduced  by  taking  a truncated  range  aperture. 

For  the  last  example  to  be  discussed,  figure  11,  a 24  Hz  infinite  CW 
source  was  located  at  a depth  of  325  ft  (99  m)  in  a deep  surface  duct  and  the 
receiver  was  deep  in  the  water  column  at  15,900  ft  (4846  m).  The  depth  of  the 
water  was  17,400  ft  (5303  m)  and  the  bottom  was  assumed  to  be  a fast  constant 
sound  speed  half  space.  The  inverse  solution  is  used  to  solve  for  the  total 
depth  dependent  Green's  function.  It  is  possible  to  take  the  solution  one  step 
further  as  was  done  in  the  previous  example  and  solve  for  the  impedance  at  the 
ocean  bottom,  but  although  the  formulas  have  been  worked  out  the  numerical  cal- 
culations have  not  been  done  to  date. 

In  discussing  the  required  range  resolution  and  aperture  (see  figure  12) 
it  is  helpful  to  recall  that  the  transfer  function  can  be  expanded  in  a normal 
mode  series.  Each  term  in  the  series  is  periodic  with  range  but  the  periods 
are  not  harmonically  related.  The  terms  which  oscillate  the  fastest  or  equiv- 
alently have  the  smallest  skip  distance  correspond  in  this  case  to  steep  angle 
bottom  bounce  paths.  The  range  resolution  must  be  chosen  small  enough  so  that 
a sufficient  number  of  spatial  samples  are  obtained  over  a complete  cycle  of 
these  paths.  For  all  of  the  results  to  be  shown  Ar  was  held  constant  at  1700 
ft  (518  m).  Because  of  the  FFT  relationship  the  extent  of  the  wave  number  ap- 
erture given  by  the  expression  in  the  upper  right-hand  corner  is  inversely 
proportional  to  the  range  resolution. 

The  maximum  range  to  which  field  values  must  be  taken  is  related  to  those 
paths  which  have  the  longest  skip  distance.  In  this  example  they  are  the  con- 
vergence zone  paths.  If  useful  information  about  these  wave  numbers  is  desired 
then  field  values  are  needed  over  a complete  cycle,  which  would  be  the  distance 
to  two  convergence  zones. 

Results  will  be  shown  for  two  different  range  apertures.  The  first  de- 
noted by  4096  points  equals  a maximum  range  in  excess  of  33  convergence  zones. 
The  second  having  256  points  extends  slightly  beyond  the  second  convergence 
zone. 


The  extent  of  the  range  aperture  then  determines  the  wave  number  resolu- 
tion as  can  be  seen  from  the  expression  in  the  lower  right-hand  corner.  The 
256  point  aperture  results  in  a larger  value  for  A?  and  thus  less  resolution. 

A comparison  is  shown  in  figure  13  between  the  absolute  magnitude  of  the 
Green's  function  for  the  exact  and  inverse  solutions  versus  the  horizontal 
component  of  the  wave  number.  For  the  range  aperture  corresponding  to  4096 
points,  the  wave  number  resolution  is  small  enough  (=  9.023420  x 10"')  so  that 
the  peaks  which  correspond  to  the  location  of  the  normal  mode  eigenvalues  are 
well  defined.  This  is  especially  true  in  that  portion  of  the  plot  correspond- 
ing to  the  bottom  bounce  (wavenumbers  less  than  approximately  .0297)  energy 
where  the  height  of  the  peaks  agree  very  well.  The  height  of  the  peaks  do  not 
agree  as  well  in  the  RSR  (wavenumbers  greater  than  approximately  .0297)  region 
but  it  is  difficult  to  see  the  difference  in  this  presentation. 

The  only  difference  between  the  result  in  figure  14  and  that  for  the  pre- 
vious figure  is  that  the  number  of  range  points  has  been  reduced  from  4096  to 
256.  Thus  the  range  aperture  has  been  reduced  from  slightly  larger  than  33 
convergence  zones  to  slightly  greater  than  2.  This  results  in  a decrease  in 
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the  wave  number  resolving  power  {*  1.443730  x lO"^)  by  a factor  of  16.  As  a 
consequence  almost  all  of  the  peaks  located  to  the  right  of  the  slide  corres- 
ponding to  RSR  paths  have  been  missed.  However,  we  are  still  doing  quite  well 
in  the  bottom  bounce  region  to  the  left.  It  is  this  region  which  interests  us 
since  the  impedance  of  the  bottom  for  these  angles  is  the  quantity  for  which 
we  eventually  wish  to  solve. 

The  previous  two  figures  provide  some  indication  of  the  error  introduced  by 
decreasing  the  range  aperture.  Additional  insight  can  be  obtained  by  comparing 
the  relative  error  between  the  exact  and  inverse  solution  at  only  those  points 
where  the  computed  inverse  solution  exist.  Those  results  are  shown  in  figure 
15  for  both  the  4096  (the  curve  which  starts  with  positive  values)  and  256 
range  apertures.  It  can  be  seen  that  the  relative  errors  associated  with  the 
two  range  apertures  are  similar  for  all  wave  numbers  and  both  do  very  well  in 
the  bottom  bounce  region. 

Thus  the  smaller  aperture  results  in  a decrease  in  the  wave  number  resolu- 
tion power  and  causes  us  to  miss  information  related  to  the  RSR  peaks,  however, 
the  solution  captures  the  peaks  in  the  bottom  bounce  region  and  also  gives 
excellent  agreement  for  those  wave  numbers  in  the  RSR  region  where  the  inverse 
solution  exist. 
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We  are  encouraged  by  the  success  of  our  results  to  date,  however,  many 
questions  remain  to  be  answered  before  the  scheme  can  be  considered  practical. 
The  effect  which  a limited  range  aperture  has  on  the  determination  of  the  im- 
pedance and  not  merely  the  depth  dependent  Green's  function  clearly  must  be 
determined.  Also  only  simulated  data  with  no  ambient  noise  has  been  used  to 
date,  and  the  real  test  of  the  practicality  of  the  inverse  solution  will  come 
when  it  is  applied  to  actual  experimental  data. 
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